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' itEtfARD AND UNCERTAINTI IN EXPLORATION PROGRAMS 
Introduction 

In our preceding paper ("Two Stochastic Models Useful in 
Petroleum Exploration")! we constructed a model to describe the results 
of wildcat drilling. Given the existence in nature of a set of targets 
of differing characteristics, in this case areal extent, we specified 
the process by which information about the targets would be accumulated. 
The model permits predictions about the success ratio and the size of 
discoveries for successive increments of drilling. Knowledge of the 
probability lavra governing the results of petroleum exploration would 
make it possible to characterize the economic risks involved, but this 
entails a more elaborate modol than the one we have proposed. It would 
be necessary, utilizing information obtained from vrells already completed, 
to specify the joint probability distribution governing the entire set 
of variables ■’hich determine the econ mic return to subsequent drilling, 

■ Previous analyses of exploratory drilling programs have 
emphasized particular aspects of uncertainty. The variable, size of 
reservoirs, has received the most attention, [Allais (1957) » Arps and 
Roberts' (1958) * Arrington (I960), Kaufman (1963).] Estimates of the 
expected value and standard deviation of reservoir size have been casually 
interpreted as measures of the economic reward and the degree of risk, 
respectively, of particular exploration programs. The size of reservoir 
found is, of course, only one aspect of the uncertainty of exploratory 
drilling. Among the other variables which have an important bearing on 
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the economics of the program are the probability of making a discovery, 
the depth of the producing formation, and. the productivity of the wells. 

In this paper vre select a set of variables wMrch are crucial 
to the economic outcome of petroleum exploration. Ttese are treated 
as random variables; the values they assume indicate; the number of 
successes that occur in a drilling program and determine, for a parti- 
cular discovery, the unit production cost and net economic return if that 
reservoir is developed. In specifying the joint probability law for 
these variables, we are forced to make extreme and probably unrealistic 
assumptions. In particular, we assume the different, random variables 
to be independently distributed, and vre do not take into account changes 
that may occur in the probability distributions as exploration proceeds. 
This latter simplification, of course, ignores the t&rust of our previous 


model, which describes a depletion process where the> largest pools, 
having been found first, are "used up" , and hence cesase to be possible 
targets, We are cpnscious of, and suitably pained. byy, s these limitations, 
and we stand ready to make .use. of- any better data-gjsnorating models 


that come along.. 

The values of the independent random variables affect the 
economic return to exploration in relatively complicated ways. As a 
consequence vre cannot deduce, the probability functions which govern 
; ; the pertinent economic measures directly from knowledge of the joint 
probability distribution of the physical variables „ Instead we rely on 
a Monte Carlo type of simulation procedure. Using, postulated probability 
functions atrl specified parameters, wo generate valines for selected 
random variables, such as reservoir size. From thl/s set of values \e 
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compute the economic magnitudes of interest, net return and unit 
production cost. This constitutes a single trial, and. the procedure 
is repeated, many times. The resulting histograms approximate the 
probability density functions of the variables which describe the 
economic outcomes of an exploratory drilling program. 

In the next section we specify the set of physical variables 
whose values are critical to the economic success or failure of an 
exploration venture. Vie then present, first, a model which relates the 
expenditures needed to develop and produce a crude oil reservoir to 
this set of variables, and, second, a model which relates the output 
of the developed reservoir to certain of the variables. Utilizing 
these models we can compute for the reservoir unit production cost 
and total value, or net economic return, the latter ■conditional upon 
the wellhead price at which the crude can be sold. 
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A Model of the Returns to Exploration 

In Table 1 wo define the variables that will be employed in 
computing the returns when a reservoir is discovered. We distinguish 
among three classes of variables! 

• (a) physical variables} which are observable upon completion 

of the wildcat well; - 

(b) certain economic variables which we postulate to be known 
with certainty; and 

(c) dependent variables} whose values we will compute. 

These variables take on a particular set of values for each wildcat 
well. 

The expenditure required to produce a reservoir is resolved 
into four components, of which the first throe comprise what is usually 
termed development investment? 

1^ " drilling investment; 

Xp ~ investment in surface gathering and processing facilities; 
Ij - camp investment, required in remote locations; and 
1^ » capitalised operating costs. 

Based on F. M. Fisher's investigation of drilling costs 

* 

[Fisher (1964)], vre assume that viell cost increases exponentially 
with depth and that; 

f* 

h = N p n (e - D + (i) 

AJ 

where 6^ is an error term. In a similar vein, it has been shown that 
the relation between investment and capacity in a chemical process plant 
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Table 1 

DEFINITIONS OF VARIABLES 


($.) physical variables 

v oil in place in reservoir, barrels 

d mean well depth, feet 

q o mean initial well productivity, barrels por day 

' . per well 


(b) economic variables 


D ratio of initial developed production capacity of 

reservoir to total proved reserves (or production 
decline rate) 

p expected price, assumed constant, dollars per barrel 

r discount (interest) rate 

T economic time horizon, years 


(c) dependent variables 

W cost of an exploratory well, dkjllars 

N number of development wells drilled 

G gross value of reservoir, present value dollars 

4 

I total expenditure required to establish and 

maintain production, develops©, nt investment 
plus operating costs, present, value dollars 

Y net value of reservoir, present value dollars 

X unit production cost, dollars per barrel 




jaay often be well approximated by the so-called "six-tenths factor 1 ' 
[Chilton ( 1930), Williams (19^7)3. Employing this we specify the 
investment for surface equipment to boj 

/J ^79 

■ h = Pa <■%> + »& N - e 2 • (2) . 

The first term on the lefthand side, with ^ 22 =0.6, represents investment 
which is dependent on the scale of operations, or throughput. The 
second component relates to expenditures which are dependent on the number 
of development wells, such as roads, gathering lines, and drilling pads, 
all of which are costly under Arctic permafrost conditions. Assuming 
that investment in a field base camp is allocated to individual pools in 
proportion to their developed capacity we have; 


% » P31 + ^3 * 


FinaOJ-y, capitalized, operating costs are assumed, to be represented by? 


A = fi 41 + V N + ^4 ■ 


Eq. (4) resolves total operating costs into a component which depends 
on the number of development wells and a component which is constant 
for the reservoir. 
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Summing tho various investment components we have: 

A-> P-i pd ' 

I ■- N [S u (e ^ - 1) 4- B 23 + 3 42 ] 

Pop 

+ P 2 1 (Nq 0 ) ^ .+ P 31 (Nq 0 ) + P 4l (5) 

+ ? 1 + ? 2 ' + ® 3 + C 4 

We assume the individual error terms comprise over a number of reservoirs 
a sequence of mutually independent random variables , identically distri« 
buted according to a known probability law. Writing u - + $ 2 + ^3 + ^4 

in Eq. (5) yields: 

I « N [P n (e X - 1) + P 23 + B 42 3 

+ P 21 (Nq o ) 22 + P 31 (Nq o ) + P 41 + U. . (6) 

The output which can be obtained from a given reservoir is 
directly related to the volume of reserves that can be proved, V , 

Postulating that the mean recovery factor, F , is known with certainty: 

* 

) 

V = F v (7) 

A number of models describing the production decline behavior of a 
reservoir have been studied, the most notable of which are exponential 
decline and hyperbolic decline. We use tho former, anticipating that 
it describes output over time sufficiently well fo-r our purposes and 
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gratefully accepting its mathematical convenience . Consequently, 
cumulative production, Q, at any given time is: 

t t 

Q t = / N t q t dt = Kq o / e " Dt dt (8) 

o - 0 

where D is the production decline rate. The term on the extreme left 
incorporates the assumption that installed capacity is not increased 
at later stages in the productive life of the pool. 

Integration of Eq. (8) yields: 

Nq n , 

Q t = -jj 2 - (1 - o'"-) . (9) 


If we consider a long period of time, such that proved reserves are 
essentially recovered V, Eq. (9) becomes: 
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Hence, for a specified decline rate, D, initial producing capacity 
can be related to reservoir size by applying Eqs, (10) and (7) I 


Nq = D F v 
o 


( 11 ) 


In addition, it can readily be seen that the number of development well 
drilled depends on reservoir size and the mean initial capacity (or 
productivity) of a well. 


V 
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N = DF 


( 12 ) 


Eq, (6) in conjunction with Eqs, (11) end (12) relates the investment 
required to produce the crude in a reservoir to the (set of physical 
parameters which describe that reservoir, as listed &n Table 1, 

We can now turn to computation of the econom&c return that 
will be gained by developing the newly discovered reservoir. With 
output as described by Eq. (8) and assuming continuous discounting, 
gross revenue must be; 




! 


0-< D «•>* dt . 


( 13 ) 


Integration yields: 


. -(rw-r)T 

o = PNq 0 C 3 


(») 


Representing the factor in brackets by A, the expression for gross 
revenue becomes; - - — -- ----- ----- — — — — 


G = P Nq o A (15) 

Net revenue, .or the eoonomic return resulting from tftie discovery of 
the reservoir, is the difference between gross revemuo and total 
investment; 

tv tv , . 

Y = G «• I (16) 
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Our basic model for the returns to exploration thus consists 
of Eqs, (16), (15)» (6), (11), and (12). This model relates the 
economic payoff of a discovery to th'e set of observable physical 
parameters in Table 1. When a wildcat is dry, the direct payoff is 
zero. The required economic parameters, also shown in Table 1, are 
assumed to be determined exogenously, to be known with certainty, and 
to be fixed for all reservoirs discovered by a series of wildcat wells. 

At the outset of exploration in a region (and therefore before the 
development of any producing capacity) the vector of parameters as 
well as the probability lavr governing u would not bo known. These 
would have to be estimated from sample data. 

Before turning to the applications of this model, it will be 
useful to describe the final dependent variable listed in Table 1, 
unit production cost, denoted by X. This will be particularly useful 
in our analysis because it is not dependent upon economic expectations; 
specifically, it is calculated without reference to expected vellhe&u 
price. By focusing on cost vo properly restrict ourselves to geological 
and technical forms of risk."*" We define unit cost, X, as the amount that 

must be realized on each barrel of crude produced in order to recover 

* 2 

the investment in the reservoir, including capitalized operating costs. 


1, M. A. Adelman (1966) has • distinguished among commercial, geological 

engineering , and political risks. Changes in selling price would represent 
commercial risk , narrowly defined. Cur estimates of expected return end 
variability of returns neglect the risk of possible changes in selling 
price and hence may be misleading. • We have not treated political risk 
here, but would do so through the revenue side, for example, by specifying 
tho probability of getting any returns after a given year. ...... 

2, This measure of cost has been used elsewhere in analyzing crude 
oil production. See Adelman (1966), Bradley (196?). 
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In the notation hb have used: 

O’ _ -JL (1?) 

- Nq 0 A 

Since the denominator of Eq. (17) represents discounted, or present 
value, output (measured in barrels), unit cost depends on the investment 
needed to obtain production and the resulting pattern of output over 

I 

| 

time. i 


1 
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ApglyxM.tho Model to the Arctic 

Iv) using this model to gain insights into the economics of 
Arctic exploration we encounter formidable problems. At the outset we 
noted the need to specify the probability laws governing petroleum 
exploration in order to properly use the data which are collected to 
make inferences about the underlying parameters. As a makeshift substi- 
tute for a moro comprehensive model we have postulated probability 
functions for tho independent random variables. In the case of reservoir 
size it is possible to employ a hypothesis which has received considerable 
testing in the literature;'' in other cases we were compelled to perform 
our own rough tests on the selected probability functions. Once the 
required probability functions are specified, we need estimates of 
parameters which characterize the Arctic area under study. Since we 
have not had access to information on which such estimates could be 

based, we conjecture possible values. It is also necessary to know the 

% 

vector of cost parameters which determine investment, denoted 8_ in the 
previous sections. To meet this need we have made some rough, calcula- 
tions which employ the estimates of another panelist, C. A, Norman, 

Tho probability functions and parameters for variables used 
in the simulations are summarized in Table 2, Parameter estimates are 
based on data describing petroleum occurrence in the Province of Alberta, 
except for initial well productivity where the Alberta data were not 
appropriate to our needs. After inspecting data for Algerian, Iran, and 

3. Two ox ample s are G. K, Kaufman (1963) and R. G. McCnossan (1969), 
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PROBABILITY FUNCTIONS AND PARAMETERS FOR INDEPENDENT RANDOM VARIABLES 
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Libya we based our rough estimates of productivity parameters on the 
Libyan information, it scarcely nsoa be said that our hypotheses 
about the distributions of the depth and productivity variables require 
further testing. It would be of great interest to test for possible 
correlation between these variables. The important figures used in 
computing the cost parameters, p, are listed in Table 3. 

With this information in hand, we return to the model of the 
previous section. The outcome of drilling a wildcat well is deter- 
mined: (a) by whether the well is a euccess, that is, whether it finds 

crude oil, and (b) if it does, by the observed values of the variables 
v, d, and q . Considering for the present only successes, we simulate 
wildcat drilling by treating v, d, and q as random variables, 

/\j i'j .... <\i . 

recognizing this distinction by the notation v, d, and q »;. A single 

A ; C.i K> 

outcome is evaluated by generating values for v, d , and q according 

i 

to the probability density functions and parameters shown in Table 2, 
•and then by computing — using Eqs. (6) , (11) » (12), (15), and . (16) — 
the corresponding values for the dependent variables listed in Table 1, 

for each set of conditions examined, 2000 outcomes were evaluated,. The 

« 

results were displayed in the form of histograms, two of which -are 

illustrated in the Appendix. Irr the next section we describe the . 
results of this procedure where the . under! ying physical and economic K-. 
parameters were chosen to represent Arctic conditions. 
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Table 3 

COKPONENT3 OF' COST PARAMETERS # 

(Derived from Estimates for the Prudhoe Bay Field) 


I Investment in producing veils: 
(a) Drilling costs, per veil 


(b) Cost increment for slant 
drilling of development, 
wells , per well 

(c) Cost of drilling pad, per well 

(d) Cost of connecting roads, 
per mile 


p 12 d 

C 1 “ B li (e " ^ 

where B-q- 119 , 000, B-j^" .0002 


(cost of drilling only, 9000 ft. 
development well, estimated to 
be <(600,000) 

^75,000 


#50,000 

#200,000 


Investment in surface facilities: 

(a) Processing plant [inc, oil and 
gas separation, gas compression 
for reinjection] 


C 2 = B a (Kq o > 


e 22 


where ^91, $ 22 “ 

(cost of plant with 100,000 bpd 
capacity estimated to bo #3.62 
million) . 


(b) Cost of injection wells , 
per well 

(c) Cost of gathering lines, 
pel' mile 

I„ Operating costs , per well per 
^ year 


I { (a,c.d) above, 

#120,000 

E^s^i+tfg/N where $ -j= 6000, 
$ 2 ~ 55^,000, and N = number of 
development wells . 


1^ Camp investment [assumed to serve #2^1,300,000 
several reservoirs in field and 
allocated to a given reservoir 
according to share of field 
output over production period.] 


^Estimates of various expenditures required to produce crude at Prudhoe Bay 
were made by C. A. Norman, We are not aware of any other estimates avail- 
able to the public which are as care full y^etailed . In adapting the original, 
figures to obtain the ones shown here andAB., ve nave combined categories and 
made simplifying assumptions for which the author of the original estimates 
should not be held responsible. 
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Simulation Results 

For the probability Functions specified in Table 2, the expected 
values for depth and veil productivity • arc , respectively, 5610 feet and 3640 
barrels per well per day. The latter figure is close to the 4000 barrels 
per well per day postulated in the Norman paper; apparently the discoveries 
at Prudhoe are deeper, however. Suppose it were known with certainty that 
ell Arctic discoveries would be at about this depth, with wells flowing at 
this output rate. The returns in such a situation- were simulated in the 
method just described, with results that are reported in Table 4. Discoveries 
under these conditions were uniformly profitable at a wellhead price of one 
dollar per barrel. 

Looking more closely at Table 4, the rows correspond to pro- 
gressively more optimistic assumptions about the size of reservoirs. At 
the low end the postulated lognormal distribution of reservoir sizes yields 
pools whose median size is around 10 million barrels (oil in place), a 
figure corresponding to experience with the attractive Devonian reefs in 
Alberta. At the high end the postulated distribution yields pools whose 
median size is about 8 times as big, a very generous assumption indeed. Pro- 
duction costs vary between 60 and 30 cents a barrel. i ho principal cause 
of cost variation lies with diseconomies of scale in surface facilities. This 
effect would be stronger, were it not that the calculations permitted a sharing 
of the costly items included under the category of base camp expenditures 
(camp, airstip, vehicles, power plant, rig mobilization) with other pools 
assumed to exist in the field. Net ro.turns increased as would be expected, 
roughly in proportion to the volume of reserves discovered. 

We now consider the situation where it is not certain that the 
development wells in the pools discovered will produce initially at 



EXPECTED RETURNS TO DEVELOPMENT 

Case I: reservoir size variable; depth and mean initial well productivity fixed 
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3600 barrels per day. Wo postulate this to be the expected rate, but 
permit ths degree of dispersion in rates which we observed in the Libyan 
data. The simulated outcomes of exploratory drilling are displayed 
in Table 5. A striking feature of this new situation is that now not 
all the discoveries are commercial --defined under our assumptions as 
being capable of producing profitably at a wellhead price of one dollar 
per barrel. The fraction of fields with positive returns is under 0.5 
for the most conservative assumption about reservoir sizes. This means 
that production cost was observed to be under one dollar per barrel 
less than half the time; the corresponding mean production cost is about 
$1.80. As would be expected, the dispersion of outcomes is much greater 
than previously; this is seen by comparing the coefficients of variation 
(defined as the standard deviation measured as a percentage of the mean). 

In Table 6 we postulate reservoir size, depth, and well 
productivity to all be variable , defined by the corresponding probability 
functions shown in Table 2, The differences between Tables 5 6 are 

not very pronounced. The coefficients of variation in Table 6 are 
generally higher; the fact that they are not universally higher suggests 
that we should increase the number of trials evaluated under each set of 
conditions beyond 2000; the observed rements of the outcome distributions 
are not yet quite stable, a consequence of the extreme skewness of the 
distributions. The means and standard deviations in the initial row of 
Table 6 are computed from the histograms shown in the Appendix:. 

The results just considered were conditional upon the exploratory 
well striking oil. To compute expected returns before drilling begins, 
these figures must be modified to take account of the probability that 



EXPECTED RETURNS TO DEVELOPMENT 

Cas© II: reservoir size, mean initial well productivity variable; depth fixed. 
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Case III: reservoir size, depth, and mean initial well productivity variable 
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the trell will in fact turn out to be a success. This might be done by 
treating successes and failures as Bernouli trials, although this clearly 
is an oversimplification. It should also be noted that the net return 
figures we have computed relate to development investment and operating 
costs, and do not include the cost of exploratory wells, which have been 
reported to run from two to three million dollars and higher on the North 
Slope, It therefore appears from Table 6 that with the figures we have 
used some combination of very favorable success ratio and large median 
reservoir size (Ugh mean of the distribution of v) is needed to make 
expected returns to a sequence of wildcats positive. Given the skewness 
of the distributions, there trill of course be some very profitable finds 
even where the expected return is low. 

• V/ith regard to the need to find largo reservoirs, for any 
existing distribution of pool sizes in nature we might expect the initial 
finds to be relatively large. This is the line of reasoning formalized by 
the model in our first paper —that the probability of finding a big pool 

i 

is higher than the probability of finding a small one. Against this, 
though, is cha benefit of better knowledge about the geology of the fegion, 
acquired as data accumulate from exploration. This might permit better 
selection among available prospects in later periods. 

4 r 

The results presented in this paper are intended to be suggestive. 
They cannot be treated ss more because we have not had Arctic data from 
which to derive our parameter estimates, and therefore have relied on 
possible similarities with already developed basins. V7e also feel that it 
will be necessary to make progress along the linos suggested in the previous 
paper before wo can confidently characterize the uncertainties of the 
exploration process. 
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Appendix 


Simulation of Met Returns 
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